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. We present a restricted solid on solid hamiltonian for fee (110) surfaces. 

It is the simplest generalization of the exactly solvable BCSOS model which 
. is able to describe a (2 x 1) missing-row reconstructed surface. We study 

£^ . this model by mapping it onto a quantum spin-1/2 chain of the Heisenberg 

q . type, with second and third neighbor SfS? couplings. The ground state phase 

ON . diagram of the spin-chain model is studied by exact diagonalization of finite 



X 



chains up to TV = 28 sites, as well as through analytical techniques. We find 
four phases in the phase diagram: two ordered phases in which the spins have 



G ■ 

O ■ a Neel-type of long range order (an unreconstructed and a missing-row recon 

O 

structed phase, in the surface language), a spin liquid phase (representing a 
rough surface), and an intermediate dimer phase which breaks translational 
invariance and has a doubly degenerate ground state, corresponding to a dis- 
ordered flat surface. The transition from the (2 x 1) reconstructed phase to 
the disordered flat phase belongs to the 2D Ising universality class. A critical 
(preroughening) line with varying exponents separates the unreconstructed 
phase from the disordered flat phase. The possible experimental signatures 

of the disordered flat phase are discussed. 
PACS Numbers: 75.10., 68.42, 64.60.C 
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I. INTRODUCTION 



Surfaces of fee metals, in particular (110) faces, display a variety of phase transitions 
which have been the subject of considerable experimental and theoretical work. 

Experimentally, (110) faces of some fee metals - like Au or Pt - are reconstructed at low 
temperature into a (2 x 1) missing-row structure, whereas other metals - like Ag, Ni, Cu, 
etc. - are not. As the temperature is raised, reconstructed surfaces tend to show two sepa- 
rate transitions: a deconstruction transition occurs first, followed, at a higher temperature, 
by a roughening transition. [|T],[| On unreconstructed surfaces, instead, only a roughening 
transition has been well characterized. |3|,|]] 

On the theoretical side, an interesting and nontrivial interplay has been anticipated be- 
tween in-plane ordering, related to reconstruction, and vertical ordering, related to rough- 
ening. H Since then many studies have been devoted to the problem. [§-|l2| 

The situation is somewhat different for the two cases, the unreconstructed and the 
missing-row reconstructed. For the unreconstructed case, den Nijs has argued that the 



phase diagram should be qualitatively similar to that of a simple cubic (100) surface. [13 



In particular, he proposes that Ag and Pd (which do not reconstruct) are good candidates 
for the realization of preroughening , an additional critical variable-exponent transition from 
a low temperature ordered phase to an intermediate temperature disordered flat phase, pre- 
viously identified in the context of restricted solid-on-solid models for simple cubic (100) 
surfaces. [|14| The unreconstructed surface possesses a Neel-type of order parameter, which 
characterizes (lxl) order 

^ixi) = ^(E^ G - r ). (i) 

r 

P(i x i) can take two opposite values at low temperature and vanishes at and above prerough- 
ening. |TT| , p!5|] (Here h r is the height of each surface atom, N is the number of (1 x 1) cells of 
the surface, and G = (2n/a x )'k is the first non-zero reciprocal vector of the lattice, see Fig. 
3.) 
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For the reconstructed case, one has four ground states to deal with (instead to the two 
ground states of an unreconstructed surface), and these can be classified, according to den 
Nijs, || by a clock variable 6 = 0, tt/2, tt, 3tt/2, see Fig. 4. The most elementary extended 
defects which one can consider are steps, which change the average height by ±1 and the 
reconstruction variable by A9 = tt/2 (clockwise or (3 x 1) steps) or A9 = —tt/2 (anti- 
clockwise or (1 x 1) steps), and Ising wall defects which do not change the average height, 
and change the reconstruction variable by AO = tt. 

The interplay between reconstruction and height degrees of freedom makes this situation 
particularly interesting. The problem has been tackled using two different lines of approach. 
The first line of approach consists in identifying heuristically the relevant defects which 
play a role into the problem, and building up models which describe their behaviour. In this 
framework, den Nijs has formulated a four state clock-step model on a length scale somewhat 
larger than the microscopic one. || On each cell of a coarse-grained lattice an integer variable 
h r - representing the average height in the cell - as well as a clock reconstruction variable 
9 r = 0, 7r/2, tt, 37r/2 are defined. By assumption, the only actors appearing in the model are 
Ising wall defects - for which (AO = tt, Ah = 0) - as well as clockwise, (AO = tt/2, Ah = ±1), 
and ant i- clockwise, (A8 = —tt/2, Ah = ±1), steps, den Nijs found that when both types 
of steps have the same energy - the so-called zero chirality limit - the model displays two 
possible scenarios: (i) if the parameters are such that the energy of an Ising wall E w is less 
than roughly twice the energy of a step E s , by increasing the temperature one goes from a 
reconstructed phase to a disordered flat phase with an Ising transition, and then to a rough 
phase with a conventional Kosterlitz-Thouless (KT) transition; (ii) in the other case - E w > 
2E S - the system undergoes a single roughening-induced deconstruction transition. || Within 
the same line of approach, but in the so-called strong chirality limit (only clockwise steps 
allowed), den Nijs |J and, independently, Balents and Kardar Jl0[ have tackled the problem 
by mapping steps into the world-lines of ID fermions, the Ising wall being represented by a 
site doubly occupied by a pair of up and down steps. The hamiltonian describing such a ID 
fermion problem is essentially a Hubbard model with an extra term describing dislocations, 
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i.e., the merging of four steps at the boundary of finite terraces on the surface. Scenario (i) 
above, realized for an attractive Hubbard U, is essentially recovered in this fermion mapping. 
As for scenario (ii), corresponding to U > 0, the roughening- induced deconstruction is now 
substituted by a doublet of successive phase transitions: first a Pokrovsky-Talapov (PT) 
transition to a "floating reconstructed" rough phase, and then a further KT transition to a 
rough disordered phase. [PlIlOl 

While these approaches are interesting per se, in that they identify potential scenarios 
for these surfaces, their outcome is obviously connected to their initial set of assumptions 
(i.e., identification of the relevant defects, and choice of their mutual interactions). 

A more direct line of approach consists in working with solid-on-solid models, which 
- while still very idealized - are fully microscopic in nature. Kohanoff et al., and, more 
recently, Mazzeo et al. have introduced a lattice solid-on-solid (SOS) model which is able 
to deal with both unreconstructed and reconstructed situations, [f^p] Their model is a 
simple modification of the exactly solvable body-centered solid-on-solid (BCSOS) model 
|TS| | obtained by adding a further neighbor interaction, and is quite closely related to the 



model we will consider here. Using a Monte Carlo simulation Mazzeo et al. study two 
separate points in its phase diagram, intended to describe the case of reconstructed Au(110) 
and that of unreconstructed Ag(110). [|ll[] For the reconstructed case they confirm a two- 
transition scenario: an Ising deconstruction transition to a disordered flat phase, followed by 
a conventional KT roughening transition. In the unreconstructed case, quite interestingly, 
they also find two transitions: a non-Ising critical one leading to a "sublattice disordered" 
flat phase, followed again by a KT roughening. [|Il[ 



The model we present in this paper is probably the simplest generalization of the BC- 
SOS model which is able to describe a missing-row (2 x 1) reconstructed phase as well the 
unreconstructed one. It essentially consists of a BCSOS model in which the sign of the, say, 
x-coupling can have either sign - the negative sign favoring the occurrence of reconstruction - 
and a further neighbor interaction in the x-direction is added in order to stabilize the recon- 
structed case against faceting into a (111) surface. The model is indeed quite similar to the 



one introduced by Mazzeo et al., except for the choice of the non nearest neighbor stabilizing 
interaction. It has the virtue of being easy to map into a one dimensional quantum spin-1/2 
model of the Heisenberg type with further neighbor SfS* interactions. Such a mapping 
provides a quite powerful way of investigating the full phase diagram of the model, without 
having to resort to Monte Carlo simulations. Moreover, the insight gained by relating our 
problem to the large body of knowledge in the field of one-dimensional quantum systems 
will appear quite evidently. 

One of the motivations of our work was a deeper understanding of the nature of the 
disordered flat phase occurring in this kind of simple modifications of the BCSOS model. 
In this respect, our study was greatly inspired by the work of den Nijs and Rommelse, in 
particular by their discussion on the many similarities between the physics of the spin-1 
and spin-1/2 chain. |14] Within our simple model, we explicitly demonstrate that, indeed, 
a spin-1/2 dimer phase is a natural candidate to representing the disordered flat phase of 
an fee (110) surface. From a surface science viewpoint, a better understanding of the phase 
diagram and of the nature of the disordered flat phases of these surfaces should help building 
future experiments directed at their study. 

Our results are summarized by the two phase diagrams of Figs. 1 and 2, pertaining to 
the spin and to the SOS model, respectively. We find four phases in the phase diagram: 
two ordered phases in which the spins have Neel-type (lxl) and (2x1) long range order 
(an unreconstructed and a missing-row reconstructed phase, in the surface language), a 
spin liquid phase (representing a rough surface), and an intermediate dimer phase which 
breaks translational invariance and has a doubly degenerate ground state (a disordered 
flat surface). The transition from the (2x1) reconstructed phase to the disordered flat 
phase belongs to the 2D Ising universality class. A critical line with varying exponents 
separates the unreconstructed phase from the disordered flat phase. Many of the above 
mentioned features are, we believe, quite robust and should hold for more complicated 
microscopic models as well. Several scenarios are possible in our phase diagram: First of 
all, the two-transition scenario for reconstructed surfaces, i.e., Ising deconstruction followed 
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by roughening; Next, the possibility - at least in principle - of a first-order transition from 
a reconstructed to an unreconstructed surface, followed by a preroughening transition to a 
disordered flat phase, and finally by roughening; Finally, for unreconstructed surfaces, either 
a single roughening transition or two transitions (preroughening followed by roughening), 
depending on the actual location of the preroughening line, which is highly model dependent. 

The rest of the paper is organized as follows. Section 2 presents the SOS model and 
its spin-chain counterpart. In section 3 we describe how the spin-chain phase diagram 
is obtained. In section 4, finally, we discuss the physics of the disordered flat phase in 
relationship to the dimer phase of the spin model, and present our conclusions. Relevant 
technical material on the spin-chain mapping is contained in the appendix. 

II. A SOLID ON SOLID MODEL AND ITS SPIN-CHAIN COUNTERPART 

The model we want to study is defined on a lattice, schematically shown in Fig. 3, which 
is comprised of two interpenetrating sublattices, conventionally referred to, hereafter, as the 
white (W) and the black (B) sublattice. In the ideal fee (110) surface the two sublattices are 
rectangular, with a x = y/2a y , and one of them sits above the other a distance a z = a y /2. The 
notation we use for the sites is such that the W and B sublattices are characterized, 
respectively, by even and odd values of i, whereas the value of j, see Fig. 3, is the same for 
a horizontal row of W sites and the row of B sites immediately above. To each site of the 
lattice we associate a height variable hij which can take only integer values (we fix a z — 1). 
Moreover, a BCSOS-type of constraint is assumed to hold, i.e., the height difference between 
each site and its nearest neighbors (belonging to the other sublattice) are forced to be ±1. 
As a consequence, h i: j are even, say, on the W sublattice and odd on the B sublattice. The 
hamiltonian for the model is written as 

Ti. = Ti-o + H-i , 

where H is the BCSOS piece 
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Ko = EEI^C^ - h wf + K 2x {h i+1J - hi.^f] , (2) 
i=i j=i 

describing next-nearest neighbor interactions with different couphng strengths in the two 
directions, and 



N X Ny 

Hi = EE^(^-M 2 . ( 3 ) 

t=i j=i 

(with K4 > 0) is a fourth neighbor interaction term in the x direction. Periodic boundary 
conditions are assumed to hold in both directions. K 2y will be always assumed to be 'positive 
and is generally the largest energy in the problem, since the y direction is the direction of the 
missing-rows we are trying to describe. The coupling in the x direction, K 2x , can instead 
have either sign as long as K4 > whenever K 2x < 0, for obvious stability reasons. A 
negative K 2x will quite clearly favour a missing-row type of reconstruction. 

The classical T = ground states for our model are easy to work out, as a function 
of the dimensionless ratio /C = K 2x /K^. One finds that for /C > the ground state is an 
unreconstructed (lxl) surface, whereas for —4 < JC < it is a (2 x 1) missing- row structure 
(with missing rows in the y-direction) , see Fig. 4. For —8 < K, < —4 a sequence of larger 
periodicity - (3 x 1 ),..., (n x 1), ... - missing-row phases is obtained, whereas for /C < — 8 
the model no longer describes a (110) surface. 

When K4 = and K 2x > 0, we recover the BCSOS model which can be exactly 



solved through a mapping to the six-vertex model, |T6[ and shows a single transition, of 
the Kosterlitz-Thouless (KT) type, between a low temperature ordered (unreconstructed) 
flat phase and a high temperature disordered rough phase. 

To proceed further, we map the model in Eqs. (0-0) into a one-dimensional (ID) quantum 
spin-1/2 hamiltonian. Details about this mapping are given in the appendix. It suffices here 
to say that the gist of the method consists in viewing the y-direction as the (imaginary) 
time direction of an appropriate ID quantum problem, whose hamiltonian H$ is selected 
in such a way that the matrix elements of the imaginary-time evolution operator e~ rHs 
coincide, in the so-called time- continuum limit, with the matrix elements of the classical 



transfer matrix of the original problem. |17| The crucial physical requirement is that the 



"time" direction must coincide with the hard direction of the classical problem, i.e., the 
coupling in the y direction has to be much stronger than the other couplings. Anisotropy is 
not expected to play a role in the phase diagram of the model, and is, moreover, physically 
quite appropriate in the present case. The nearest-neighbor height constraint suggests that 
the relevant degrees of freedom can be described by introducing a spin- 1/2 Hilbert space, 
and mapping the height difference between nearest-neighbors within the same "time-slice" 
- the dashed line in Fig. 3 - to the z-component of the usual spin operator 

S! < — > ^(hi+1,0 - hi i0 ) . (4) 

Fig. 4 illustrates explicitly the mapping of the (lxl) and (2x1) ground states in terms of 
spin configurations along a time-slice. In terms of spin operators, the quantum hamiltonian 
reads 

J Nx Nx J 2 

Hg = — — 22[Si~S i _i + Si SiLi] + 22[JzS*Si-i + J2S-S-_ 2 + TrS'Sts] > (5) 
A i=i i=i / 

where the spin couplings are related to the original couplings as follows: 

J = 2exp{-4(3K 2y ) 
J z = 8P(K 2x + 3K 4 ) 

J 2 = 16(3K 4 , (6) 

and f3 = l/fe^T. Periodic boundary conditions on the heights are inherited by the spin 
model, and also imply that we need to work in the spin sector with zero total magnetization. 
It is well known that the mapping is such that the free energy per site of the classical problem 
is related to the ground state energy per site of the one-dimensional quantum problem, i.e., 



(3f = €gs- |0§ The temperature clearly enters through the spin couplings, see Eqs. |6|, 
so that any temperature singularity of the classical free energy can be seen as a ground 
state energy singularity for the quantum problem as a function of the couplings J z j J and 
J2/J. Moreover temperature averages for correlation functions of the classical problem can 
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be likewise rewritten in the form of ground state averages for the corresponding quantum 
correlation function. QTSJ In summary, to obtain information about the temperature phase 
diagram of the classical model we need to study the ground state phase diagram of the 
spin-chain model. 

As it is, the spin model in Eq. |5] has not been previously studied, to our knowledge. 
A similar model, however, has been discussed in Ref. ||19|| , for instance, and has been the 
topic of a quite extensive literature. It differs from our model in that the second neighbor 
coupling is taken to be spin isotropic, and the third neighbor coupling is missing, i.e., 

j N x N x 

H = -i Yltfsu + s-sti) + T\J z stsu + j 2 s 1 ■ sU] . (7) 
z i=i i=i 

We will see later on that many features of the phase diagram of the model in Eq. [7| are 
indeed present in the phase diagram of our model too. The only exception is an extra phase 
appearing in our model for sufficiently large J2, a phase which corresponds, in the surface 
language, to an ordered (2x1) missing-row structure. 



III. ANALYSIS OF THE QUANTUM SPIN-CHAIN MODEL 

We start this section by defining the interesting order parameters and correlation func- 
tions that we need to look at in order to determine the phase diagram for our model. A quite 
direct way of defining an order parameter for the quantum spin-chain case, is to consider the 
"square" of the appropriate "staggered magnetization" A -1 J2j e % ^ p ^Sj. The result can be 
trivially cast in terms of the Fourier transform S zz (k) of the spin-spin correlation function 
(the "structure factor"): 

J?xi = Hm 1 £ e^(S^) = Jim ^(tt/p) , (8) 
3 

where the average has to be intended as a ground state average. Clearly, ordered phases will 
be signalled by a structure factor Sff(k) which diverges linearly with N - for the appropriate 
k-vector - as the length of the system N — > 00. More specifically, S^(tt) will diverge for 
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a phase with |||| (Neel) long range order (LRO), a (1 x 1) unreconstructed phase in the 
surface language, see Fig. 4. A divergence of S^(n/2) will instead indicate |tl! LRO, i.e., 
a missing- row (2x1) reconstructed surface, see Fig. 4. 

Another quantity one needs to look at, in a SOS model, is the height-height correlation 
function 

G(r - r') = {[hr - h T >f) , (9) 

which distinguishes a flat [G(r) < const as r — > oo] from a roug/i [G(r) ~ In (r) as r — > oo] 
phase. The height-height correlation function has a simple translation in the spin-chain 
model. Restricting our consideration to correlations within the same horizontal row, say 
j=0, we have 

rt— 1 n 

G(n) = ([h nfi - h 0fi ] 2 ) = 4 £ (GS\S?S!\GS) = n + 8 J> - i){GS\S^\GS) , (10) 

ij'=0 i=l 

where the last equality follows from assuming translational invariance for the spin-spin 
correlation function. Notice, in passing, that the term linear in n is exactly cancelled by 
working at zero total magnetization since, in that case, 

52(GS\SSS?\GS) = -(GS\S Z S Z \GS) = -- . 

Moreover, it is easily checked that whenever the spin-spin correlation function (SqS*) pos- 
sesses a large distance uniform term of the type —K/(2n 2 n 2 ), G(n) will diverge logarithmi- 
cally as 

AK 

G{n) = — log (n) + • • • n^oo, (11) 

signalling a rough phase. 

We now discuss the information that one can get on different regions of the phase diagram 
from analytical approaches. In the limit J z /J — > oo and J 2 / J — > oo things are very easy 
to work out. Setting J = 0, one finds that for J 2 < (2/3) J z the Neel states | fltl • • •) 
and | | TIT ' " ') are the two possible ground states, whereas for J 2 > (2/3) J z there are four 
possible ground states 
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(12) 



Such a fourfold degeneracy of the ground state reflects the fourfold degeneracy of the (2x1) 
missing-row surface, as illustrated in Fig. 4. 

We now turn to the region of the phase diagram close to J z = J2 = 0. Clearly, when 
J2 = (i.e., K4 = 0) the problem reduces to the XXZ Heisenberg chain - the quantum 
spin-chain counterpart of the BCSOS model - whose physics is well known. For \J Z /J\ < 1 
the model is gapless and critical, i.e., the ground-state spin-spin correlation functions decay 



algebraically at large distances p9|,[20[] 



(5 *5*>«cos(7rn)A__A_ + ... 

(5 + 5->« cos (rrn)^ + ■■• . (13) 

The exponents rj and r\ z are exactly known: 

77 = — = — arccos {—J z l J) ■ (14) 

From the previous discussion - see Eq. O- it follows immediately that such a phase represents 
a rough surface. For J z j J > 1 the model has a gap and a Neel order parameter, and 
represents the already mentioned unreconstructed (1 x 1) phase (the BCSOS flat phase). 
The transition at J z /J = 1 is of the KT-type (the BCSOS roughening transition). Consider 
now the case J2 > 0. By applying a Wigner- Jordan transformation we can equivalently 
rewrite our spin model in terms of a spinless fermion model 

BZ 1 BZ 

H f = -JY: cos (k) c+c k + +-£ V(q)p(q)p(-q) , (15) 

k 1 

where p(q) = J2k c t c k+q is the density operator and V(q) = J z cos(q) + J2COs(2g) + 
( J2/2) cos (3g) is the interaction potential. Zero total magnetization for the spin model 
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implies a half-filled band for the fermion model. Applying standard techniques in the theory 
of ID fermionic systems, i.e., linearizing the band around the two Fermi points at kp = ±7r/2 
and keeping only interaction processes around the Fermi points, we can write down an effec- 
tive continuum field theory ("g-ology model"). [^2j Let a± ) k be the destruction operator for 
right-moving (+) and left-moving (— ) fermions around the two Fermi points, and let p±(q) 
be the corresponding density operators. The g-ology model reads: 

H e s = v F Y^ M a +,fc a +,fc _ a -,k a -A 

k 

+ T^P+( q ">P+(~ q "> +P-(<?)P- (-<?)] + 



+ y Eb (+) (<?)p (+) (-9) + p h (?)p h (-?)] 

L q 



(16) 



where p^ + \q) = Y,k a +,k- q a -,k an d P^~'{q) = Y.k a ~.k- q a +,k- To lowest order, the couplings 
turn out to be g^ = J z + (3/2) J 2 , g 2 = 4J 2 + 2 J 2 and g% = —2J Z + J 2 . The g 3 term describes 
umklapp processes. In absence of g 3 , the g 4 and g 2 terms define a spinless fermion Luttinger 
model which can be exactly solved (at least when the band cut-off is removed and one works 
with infinite linear bands). Notice that J 2 competes with J z in determining the sign of 
the umklapp term coupling constant g%. This is a quite general condition under which a 
preroughening line, separating two different gapped phases, is obtained. The continuum 
model coincides, apart from slightly different weak coupling values for g 2 ,(?4, and g%, with 
the corresponding continuum model derived for the hamiltonian with spin isotropic second 
neighbor interaction in Eq. [?]. Jl9| The physics of the two models should indeed be the same 
in the region close to the origin, and we can therefore borrow most of the conclusions of 
Ref. [19]. Summarizing the results of Ref. [l]|, we can say that there is a whole region of 



the phase diagram - for sufficiently small J 2 and J z - in which the model is within the basin 
of attraction of a spinless Luttinger model, i.e., g% scales to zero upon renormalization. In 



short, borrowing Haldane's terminology, [23] we have a Luttinger liquid. As a result, the 



spin-spin correlation functions will have exactly the same large distance behaviour as in 
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Eq. [13], with suitable coupling dependent exponents r\ z and rj. Moreover, quite generally, 
r] z = I/77 and the constant K is related to the exponents, K = rj z /2. [^] In the presence 
of the umklapp term (#3 7^ 0), such a Luttinger liquid becomes unstable towards phases 
having a gap in the excitation spectrum whenever the constant K goes to 1/2 (or rj z — > 1). 
2~2"| Such a transition is predicted to be in the KT universality class. There is nevertheless 



a whole line of points, inside the gapped phases, on which effectively the coefficient of the 
umklapp term g% is exactly zero. Such a line will consist of critical points in which the model 



still behaves as a Luttinger liquid with a varying exponent r\ z < 1. The theory of Ref. [19 
predicts also the nature of the two gapped phases to which the Luttinger liquid becomes 
unstable on either side of the " g 3 = 0" line. On one side, the already mentioned Tit! Neel 



phase is found, whereas on the other side a spontaneously dimerized phase should occur. [19 



To determine the actual location of the different lines in the phase diagram we need 
to resort to a numerical study. We have performed a finite size scaling analysis of exact 
diagonalization data for chains up to N = 28 sites. The location of the KT lines terminating 
the Luttinger liquid (rough) phase are determined from the size scaling of the gaps as done, 
for instance, in Ref. |2U|| . More in detail, the finite size gaps from the ground state to the 
first excited state in the spin-1 and spin-0 sectors will determine r\ and rj z , respectively, 

N[E N (S=1, k = 0) - E N (S=0,k=0)) = 7rv s r] + ■ ■ ■ 

N[E N (S = 0,k = n) - E N (S=0,k = 0)] = Trv s r] z + --- , (17) 

as iV — > 00. The sound velocity v s can be similarly determined by knowing the gap from 
the ground state to the first excited state of smallest momentum 

N[E N (S = 0, k = 2n/N) - E N (S = 0,k = 0)} = 2ttv s + ■■■ . (18) 

An independent way of estimating rj z (or K) comes from the q — ► limit of the structure 
factor Sjf(q). In the Luttinger liquid phase we have, for q — > 0, S zz (q) = (K/2n)q + o(q), 
and therefore 

lim NS z N z (q = 2n/N)=K . (19) 

N^oo 
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(The advantage of the latter relationship is that the sound velocity v s does not appear 
at all.) As an illustration, Fig. 5 shows the finite-size results obtained for v s tj (a) and v s 
(b) along the line at J z /J = 1 for several values of J 2 /J- Power law extrapolations of 
data for v s r], v s r] z , and v s yield the results for the exponents summarized in Fig. 6. Quite 
unambiguously, a transition - at which i] z = 1/rj = 1 - occurs for J2/J ~ 1-2 -f- 1.3. An 
independent confirmation of such a result comes from the structure factor data shown in Fig. 
7 (see Eq. |19|), from which it is quite clear, for instance, that J2/J = 1.4 is already beyond 
the transition point (the dashed line in Fig. 7, at K = i] z /2 = 1/2). A similar procedure 
has allowed us to determine the location of the KT lines in the phase diagram of Fig. 1. 

The location of conventional second order transitions (which are expected to terminate 
the large J 2 ordered phase) are determined from the "specific heat" C v = d 2 eGs/dJ 2 2 . Fig. 8 
shows C v for chains up to iV = 28 sites as a function of J 2 at J z / J = 1. The location of the 

(c) 

maximum of C v quickly saturates to our estimate of the critical value of J 2 , J 2 ~ 1-66 J. 
Plotting C v at the estimated versus the logarithm of the system size (see inset of Fig. 8), 
we conclude that C v diverges logarithmically, i.e., the critical exponent a = 0. In order to 
calculate another exponent, we consider the structure factor S z ^{n/2). In the ordered (2 x 1) 
phase Sn(7t/2) diverges linearly with N. Exactly at the critical point, however, the spin-spin 
correlation function will decay as a power law at large distances, (SqS^) ~ cos {rm /2)/n v . 
If the power law exponent rf happens to be strictly less than 1, the structure factor will 
diverge at the critical point as 

S z n z (tt/2) w (const) N 1 ^' at J 2 (c) as iV -> 00. (20) 

Fig. 9 shows a log-log plot of S^(tt/2) versus iV for J z = J and different values of J 2 around 
~ 1.66 J. A good straight line fit is obtained only for J 2 = 1.66 J, and the slope of 
the straight line is ~ 0.75 (the solid line in Fig. 9). This implies rf = 1/4. The transition 
belongs, therefore, quite unambiguously to the 2D Ising universality class. 

As far as the predicted dimer phase is concerned, there are clear numerical indications 
that we are indeed dealing with a doubly degenerate ground state with a gap above it. Fig. 
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10 shows, for instance, the size dependence of the gap to the first excited S = 1 state for 
J2/J = 1.6 and J z /J — 1. A pronounced increase is seen when one plots NAE^ versus 
1/N, a quite clear hint that AEn is saturating to a constant. On the contrary, the gap to 
the first excited S = state of momentum 7r is closing up, as shown in the inset of Fig. 10, 
hinting that such a state is actually degenerate with the ground state in the infinite volume 
limit. Notice that, with chains of length iV up to 28, a reasonably clear gap, as the one 
illustrated in Fig. 10 for the S = 1 state, is seen only quite far from the transition to the 
dimer phase. The same argument applies to a direct test of LRO. We have in fact measured 
a dimer-dimer correlation function 

D ij = ((SUSfKS^S*)) , (21) 

which is predicted to acquire LRO, m C(— l) i— J , in the dimer phase. Once again, 
although a clear tendency to LRO is seen within our "small chains" , a full fledged signal of 
LRO, i.e., the k = it structure factor of diverging linearly with N, is almost never seen, 
except very far from the transition. All these features point into the direction of small gaps 
and large correlation lengths, which are numerically difficult to assess. 

The location of the variable exponent (#3 = 0) line in between the Neel and the dimer 
phase - characterized by 1/4 < i] z < 1 |19[ - has been obtained from a combination of 



the gap scaling procedure in Eq. [17], and from the power law divergence of the structure 
factor Sjv (71") (using the procedure discussed above for assessing the Ising nature of the 
second order transition). The variable exponent line, labelled PM in Figs. 1 and 2, can be 
identified with the preroughening line introduced by den Nijs and Rommelse in the context 
of restricted SOS models for simple cubic (100) surfaces. [[L4[] The preroughening line merges 



with the second order (Ising) line at a tricritical point M, beyond which a first order line 
should follow. Such a first order line should asymptotically tend to approach J 2 = (2/3) J z , 
the value determined from the J = limit. 

The final phase diagram for the spin model is depicted in Fig. 1. Most of the lines 
are quantitative although some points are not determined with great accuracy. The region 
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around the point labelled P shows pronounced size effects which make the exact location 
of the lines somewhat more uncertain than elsewhere. Even more pronounced size effects 
are seen is the J z ~ region of the phase diagram where the KT line ending the spin-fluid 
(Luttinger) phase approaches the Ising deconstruction line. Such large size effects, if on the 
one hand prevent us from giving a clear answer to the behaviour of the model in that region, 
are, on the other hand, a quite probable hint that the Ising and KT lines might merge at 
some point. We are at present unable to say if the scenario of a single "roughening-induced 
deconstruction" is indeed at play in our model, and even less able to say if, in such a case, 
the Ising and KT nature of the original lines would simply be superimposed, as seen in Ref. 

A translation of the spin-chain phase diagram in Fig. 1, using the coupling relationships 
in Eq. |6|, allows us to determine the temperature phase diagram for the initial classical SOS 
model, which is shown in Fig. 2 for K^jK% y = 0.1. The different lines can be obtained, 
in principle, in a quantitative way from Fig. 1, with discrepancies from the expected ex- 
act results which are probably around few percents, even for not so strongly anisotropic 
couplings. 



IV. THE DIMER PHASE AS A DISORDERED FLAT PHASE 

Let us summarize some well known properties of the dimer phase appearing in our phase 
diagram in Fig. 1, trying to give an interpretation of them in the surface language. (See 



also Ref. [fL4[ .) The dimer phase has a doubly degenerate ground state with a gap above 
it. Translational invariance is spontaneously broken. All two-spin correlation functions 
decay exponentially to zero at large distances, which implies - see discussion in section 2 - 
that the surface is disordered but flat. There is, however, LRO if one considers four-spin 
correlation function of the type ((Si-i ■ Si)(Sj-i ■ Sj)). To be more specific, let us consider a 
representative dimer phase point which corresponds to the exactly solvable case of the model 
in Eq. [7], i.e., J z = J and J 2 = J/2. This is a particular case of a large class of spin-chain 
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models with competing interactions whose ground state is known to be a simple product of 
singlets (Madjumar-Gosh models). |2^] Explicitly, for any finite (even) size N there are two 
ground states (which go into one another by translation of a lattice spacing) 

= |12)|34).-.|JV-1JV) 
|* 2 ) = |23)|45) ■■■\N- 2N- l)\Nl) , (22) 

where \ij) = | fl — IT)/v2 denotes a singlet between sites i and j. Expanding the product 
of singlets it is easy to check that each time a pair f j appears, it must be followed, sooner 
or later, by a J, j pair. The surface is therefore flat since a step up is followed necessarily by 
a step down. 

Let us now ask ourselves the following question. Suppose we want to count, in the 
surface terminology, the difference in the number of white and black local maxima in the 
surface. For simplicity, we restrict our considerations to sites which are local maxima when 
considered in the x-direction only. In the spin language, a local "maximum" at site j 
occurs whenever the site j — 1 has spin j and the site j has spin j. An operator which 
"counts" the maximum at j is therefore {Sj_ l + l/2)(l/2 — Sf). The difference between 
white (even j) and black (odd j) maxima is therefore measured by the order parameter 
P BW = {2/N)J2 j e i7Tj {S^ l + l/2)(l/2 - SJ). P BW is odd under translation. Its value 
is 1 on the Neel state | T-lt-l "'•)> an d —1 on the other Neel state | HIT ' ' ')• (Quite 
generally, it is different from zero in the whole Neel phase of our phase diagram in Fig. 1.) 
Consider now the value of Pbw 011 the dimer state Using the elementary results that 

(^ r i|S'|-_ 1 S'|-|^ r i) = —1/4 and (^ r i|S'|-S'|- +1 |^ r i) = 0, and observing that the terms linear in 
S z vanish identically, we arrive at 

(*i|^W|*i> = ~^E^ j (nSU S j\^) = l ■ ( 23 ) 

(Similarly, (^I-PbuH 1 ^) = —1/4.) We clearly recognize, in Eq. ^3], the zz-component of 
the usual dimer order parameter. Therefore, a surface interpretation of the dimer order 
parameter and of the spontaneous breaking of translational symmetry is that either the 
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white or the black sublattice tend to dominate in the surface local maxima. This property 
has been directly verified by a Monte Carlo simulation of the classical SOS model. |2IJ More 
precisely, we have considered an operator Oij whose value is 1 if the site is an actual 2D 



local maximum, and zero otherwise, and then constructed the following order parameter [25 



P ( bw } = T^ r (E(-l) i OM) • (24) 

Our Monte Carlo simulation confirms that Pr^ 7^ (in the thermodynamical limit) in the 
disordered flat phase of Fig. 2. 

V. ORDER PARAMETERS AND SURFACE SCATTERING EXPERIMENTS 

As exemplified in the previous section, a disordered flat phase corresponding to the 
spin dimer phase should be characterized by a spontaneous breaking of the translational 
invariance, and, consequently, by a predominance of one of the two sublattices in the topmost 
layer. 

These features should be of relevance in the context of surface scattering experiments. 
We discuss here the case of He scattering. The discussion of the X-ray scattering case can 



be carried along similar lines, and is dealt with in Ref. p5" |. In the kinematical approxima- 
tion, and within a SOS framework, the coherent part of the specular peak intensity with 
perpendicular momentum transfer in the so-called anti-phase configuration is proportional 
to 

r oh (Q = 0,q z = ir/a z ) oc | (£ e^-c^) | 2 , (25) 

where oiij is an appropriate shadowing factor which takes into account the physical re- 
quirement that peaks scatter more then valleys. [11,26] For our BCSOS-type of model, in 



which even z's (W sublattice) are associated to even heights /iy and odd z's (B sublattice) 
to odd hi j, one immediately concludes that e l7Thl < j = (—1)* for any allowed height config- 
uration. I coh (Q = 0, q z = 7f/a z ) would therefore be identically zero if all the atoms were 
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to scatter in the same way (ay = 1 for all In the assumption that only the local 

peaks scatter efficiently (ay = 1 if is a local peak, ay = otherwise) we obtain that 
I coh (Q = 0, q z = ir/a z ) is proportional to the square of the order parameter introduced 
above, 

I coh {Q = 0,q z = ir/a z ) oc \pffi\* . (26) 

Quite generally, for a reasonably large class of choices of shadowing factors ay, the breaking 
of translational invariance should guarantee that I coh (Q = 0, q z = ir/a z ) is different from 
zero (albeit possibly small) in the disordered flat phase considered here. More precisely, this 
is so for all the shadowing factors which can be written in terms of local operators (in the fry 
variables) whose correlation function is long-ranged in the disordered flat phase. We mention 
here a particularly simple choice of shadowing factors, proposed in Ref. ||26|| , which does not 



involve long-ranged operators: ay = 2 — riy/2, where ny is the number of neighbors of 
the atom in which are found at a level higher than the atom itself. This expression 
linearly interpolates between a = 2 (local maximum) and a = (local minimum). It can be 
recast in the form ay = 1 — (1/4) Xm.n. A/iy, where — A/ty is the height difference between 
site and any of the four neighboring sites. Indeed, by exploiting this linearity, it is 
very simple to show that, such a choice of a»j leads to a I coh (Q = 0,q z = ir/a z ) which is 
proportional to the square of the (1 x 1) order parameter Pn x i) (see Eq. [I]), 

I coh (Q = 0,q z = 7r/a z ) ex |<£(-1)7^.)| 2 , (27) 

and therefore vanishes at and beyond the preroughening line. 

In conclusion we have proposed a simple microscopic SOS model for fee (110) surfaces 
which can be easily mapped into a spin-1/2 quantum chain of the Heisenberg type with 
competing SfSj interactions. The resulting phase diagram is quite simple and clear. Apart 
from the obvious unreconstructed and missing-row reconstructed ordered phases, and the 
high temperature rough phase, we find a disordered flat phase which shows the physics of the 
dimer spin phase. The possible experimental relevance of such a state has been discussed. 
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APPENDIX: DERIVATION OF THE SPIN MODEL 



The general idea behind the mapping of a classical statistical mechanics problem in D 
dimension into a quantum problem in D — 1 (space) dimensions is quite old, and needs not 
to be repeated here. |T3 For the problem at hand, the details of the mapping are quite 



similar to the procedure sketched in Ref. [13 for the spin-1 case. We briefly report them 
here for the reader's convenience. 

The starting point for the mapping to a spin model is a T-matrix formulation for the 
classical partition function 

Z = Y J t- m = Y.- ■ ■ E (h (1) \f\h^)---(h^\f\h^) ,(h®\f\hW) , (28) 

hi,j hW h^y) 

where \h>") = {hij : i = 1, • • • , N x } is the jth row configuration (the dashed line in Fig. 3), 
and T is the classical transfer matrix. Periodic boundary conditions have been assumed. It 
is implicitly assumed that the configurations included in the sum have to respect the BCSOS 
restriction constraint. Indeed, by exploiting this feature, which implies that hi+xj — hij = 
±1, we can associate to any row configuration \h^) a state \j) = \Sf , Sf , ■ • ■ , Sfj x ) in the 
Hilbert space of a quantum spin-1/2 chain (of length N x ) by the relationship 

Si <— > \{hi+ij - tkj) . (29) 

(In so doing we actually loose information on the absolute height of the surface.) The idea 
is now to try to reproduce the Boltzmann factors appearing in the matrix elements of the 
classical transfer matrix (h^ +1 '\T\h^>) by a suitable quantum operator Tq in the spin Hilbert 
space, i.e. 

(hV +1) \f\hP) = <y + l\T Q \j), (30) 

where \ j) and \j + 1) are the spin states corresponding to \h^') and \h^ +1 ') respectively. For 
the hamiltonian in Eqs. (0-0), the T-matrix element reads: 

N x 



(hV+Q\f\h®) =exp{-YlK 2x (h% 1 -h^rf + K i (h% 2 -h^ 2 )*]} x 

i=i 
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N x 

exp{-^E(^' +1) -#) 2 }- (31) 
i=i 

The first exponential contains interaction terms within row j and it is quite simple to rewrite 
it in terms of exponentials of S z operators only (see below). On the contrary, the second 
exponential contains the interaction which couples row j to row j + 1 and is explicitly off- 
diagonal. It is not hard to show that, given all the height contraints, the spin operator 
{SfS^ + S~Sf_^) does exactly the job of transforming the spin configuration around site i 
in row j, with h!f \ into the correct spin configuration of the corresponding site in row j + 1 
with h}j +V> = h^p ± 2. Such an operator must therefore be multiplied by e~ 4K2 y if we want 
to reproduce the correct Boltzmann weight. The (diagonal) situation in which h^ +1 ^ = 
is simply given by the identity operator in the spin space. The only tricky point is that 
(S' i f S'. i l 1 + SfSiLi) does not commute with the pieces involving S z , so that care must be 
taken in ordering the various terms. One can show that the operator Tq which does the job 
correctly is given by 

T Q = ■ ■ ■ T{n + 3)T(n) • • • T(7)T(4)T(1) 
T{n) = [1 + e-V^iS+S^ + 5"5+ J] x e -^^ n _ 1+m 

Xe -4/^4(s*_ 2 +s*_ 1+ s*+s* +1 ) 2 ^ ( 32 ) 

where every T(n) for n — 1, • • • , N x is included in the product which is however organized 
in a cyclic way (i.e., with n increasing by 3 from right to left and taken modulo N x ). This 
is not yet a very useful relationship, because, as mentioned above, the quantum operators 
appearing in the above expression do not commute. Our aim is indeed to work in a limit 
in which all the different factors can be grouped together under a single exponential of the 
form e~ Hs . This is the case in the so-called time continuum limit, defined by the fact that 
e -WK 2y = is sma ] 1> in which case i + (J/2)(. • •) w e ( J / 2 )(-), and \(3K 2x \ as well as \(3K A \ 
are also small, so that one can neglect the commutators coming about in regrouping all the 
exponentials. The final result can be cast in the form 

f < — > T Q w e~ Hs (Time Continuum Limit) , (33) 
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where the spin hamiltonian H$ is given in Eq. [5|. The time continuum limit is essentially a 
limit of high anisotropy. Whenever anisotropy is not "relevant" , one obtains a temperature 
phase diagram for the classical problem from studying the ground state phase diagram of the 
quantum problem. What is more surprising is that the temperature phase diagram turns 
out to be quantitatively accurate (within a few percent) even for situations in which the 
classical system is not so strongly anisotropic. This can be explicitly checked in the 2D Ising 
18H and in the BCSOS case. 
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Figure Captions 



Figure 1 Ground state phase diagram for the spin-chain model. PM is a variable exponent 
(preroughening) line with 1 > r\ z > 1/4. Beyond the tricritical point M the character 
of the line should change to first order. Ground state degeneracies are given in square 
brackets. The dashed line indicates the region where strong size effects prevent us 
from getting accurate data for the exponents. 

Figure 2 Ground state phase diagram for the SOS model as obtained from Fig. 1 using the 
coupling relationships between the two models. We have taken here = 0.1K 2y . 

Figure 3 Schematic top view of the lattice. The site notation is explicitly indicated. 
The white and black sublattices are denoted by open and solid circles, respectively. 
The dashed line represents the time-slice set up for performing the spin-chain mapping. 
The dimensions (a x , a y ) of the (1 x 1) unit cell are shown. 

Figure 4 Schematic height profiles, along a time-slice of Fig. 3, of the two ground states 
of the unreconstructed surface, and of the four missing-row ground states with their 
clock label. The corresponding spin configurations are explicitly indicated. 

Figure 5 (a) Size scaling of the gap to the first S = 1 excited state, (b) Size scaling of the 
gap to the first excited state of smallest momentum. 

Figure 6 Summary of the results for the exponent rj and r/ z along the J z = J line. A 
transition, at which rj = l/rj z — > 1, is clearly seen for J2/J ~ 1.2 -j- 1.3. 

Figure 7 Size scaling of the q — > limit of the structure factor for several values of J 2 along 
the J z = J line. 
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Figure 8 Size scaling of the "specific heat" C v = d 2 tGs/dJ% along the J z = J line. The 
peak position saturates at ~ 1.66J. The logarithmic nature (a = 0) of the 
divergence of C v is demonstrated by the inset. 

Figure 9 Size scaling of the %/2 component of the spin structure factor showing the iV 3//4 
divergence at the Ising critical point ~ 1.66 J. 

Figure 10 Size scaling of the gap to the first excited S = 1 state in the dimer phase 
(J 2 /J = l.Q,J z /J = 1). The inset shows that the S = state of momentum tc is 
degenerate with the ground state for iV — > oo. 
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